Potential energy curves, transition dipole moments, and radiative lifetimes of KBe molecule
Wan Ming-Jie1, †, Jin Cheng-Guo1, Yu You2, Huang Duo-Hui1, 3, Shao Ju-Xiang1, 3
Computational Physics Key Laboratory of Sichuan Province, Yibin University, Yibin 644007, China
College of Optoelectronic Technology, Chengdu University of Information Technology, Chengdu 610225, China
School of Physics and Electronic Engineering, Yibin University, Yibin 644007, China

 

† Corresponding author. E-mail: wanmingjie1983@sina.com

Abstract

An ab initio calculations on the ground and low-lying excited states (X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+) of KBe molecule have been performed using multireference configuration interaction (MRCI) plus Davidson corrections (MRCI+Q) approach with all electron basis set aug-cc-pCV5Z-DK for Be and def2-AQZVPP-JKFI for K. The 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+ states are investigated for the first time. Inner shell electron correlations are computed on the potential energy curves (PECs) calculations. The spectroscopic and molecular parameters are also predicted. In addition, The transition properties including transition dipole moment, Franck–Condon factors qυ′υ″, Einstein coefficients Aυ′υ″, and the radiative lifetimes τυ′ for the , , and transitions are predicted at the same time.

1. Introduction

In the last twenty years, a number of interesting concepts have been discussed in cold and ultracold molecules formed from alkali-metal atoms (Alk) of group I and alkaline-earth-metal atoms (AEM) of group II. Photoassociation (PA)[1] is one of the most common approaches to obtain ultracold molecules from ultracold atoms. PA spectroscopy of KRb,[2] Yb2,[3,4] YbRb,[5] and RbCs[6,7] has been observed.

Spectroscopic parameters, rovibrational levels, transition dipole moments, and Franck–Condon factors are also important data in the study of interactions between ultracold molecules.[8,9] Recently, several theoretical investigations have focused on the electronic structures of ultracold molecules,[822] such as MgAlk molecules,[9] SrAlk molecules,[10] AEMLi molecules,[11] KRb,[12] LiBe,[13,14] NaBe,[13,15] KLi,[16] KBe,[17,18] RbBe,[19] CsBe,[19] RbYb,[20] LiYb,[21,22] and our previous works for MgCl,[23] MgBr,[23] BeCl,[24] BeBr,[24] BeI,[25] MgI.[25] There are a few theoretical investigations on the potential energy curves (PECs) of KBe molecule. Bruna et al.[17] have given the equilibrium distance and dipole moment for the ground state of KBe molecule for the first time. The PECs of four doublet states with different inner electron correlations were investigated by Xiao et al.[18] Four sets of the frozen core orbitals (FCOs) for each electronic state were considered in their work. At first, they kept the 1s orbital of Be and 1s2s2p3s3p orbitals of K as frozen core orbital; then, 1s orbital of Be, 1s2s2p3s orbitals of K; 1s orbital of Be, 1s2s2p orbitals of K, and 1s2s2p orbitals of K were considered FCO. They believed that the spectroscopic parameters based on the last one were more reliable.

In this work, we have carried out the multi-reference configuration interaction (MRCI)[26,27] computations on three , two 4Π, and three states of KBe molecule based on inner shell electron correlations. At the present time, Li,[28] Zhu,[29] Wang,[30] Liu,[31] and Xing[32] have adopted the MRCI method to calculate the ground and low-lying states of CS+, AsO+, P2, GeO, and CF+, and the obtained results conform to the experimental values. No experimental data have yet been found in the literature for these states, and our main goal is to provide accurate spectroscopic parameters for these states. In addition, the dipole moments, the transition dipole moments (TDMs), Einstein coefficients (Aυ′υ″), and Franck–Condon factors (qυ′υ″) are also presented.

2. Computational details

The electronic structure calculations are performed with the MOLPRO 2010 program package.[33] Working in C2ν point group symmetry, which has four irreducible representations (A1, B1, B2, and A2). The A1 irreducible representation yields Σ+ states and a component of Δ states; B1 gives the Π states, and A2 yields Σ states and the other component of Δ states.

Complete active space self-consistent-field (CASSCF) calculations are carried out. Using the CASSCF[34,35] wave functions as a zero-order function, the energies of the Λ–S states are computed via the MRCI method plus Davidson corrections (MRCI+Q).[36,37] Eight electronic states (X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, 34Σ+) are investigated in the present work. The PECs for the , , and 4Π states are calculated separately. To improve the quality of PECs, core–valence correlations are computed in PECs calculations in MRCI+Q step. In CASSCF step, the 1s2s2p orbitals of K are doubly closed, which are frozen orbitals, 1s of Be and 3s3p of K are closed-shell orbitals, the 2s2p orbitals of Be and 4s4p3d orbitals of K are the active orbitals. Considering electron correlation effect in the MRCI+Q step, the 1s2s2p orbitals of K are core orbitals, 1s of Be and 3s3p of K are closed-shell orbitals, therefore, 13 electrons are correlation electrons.

In the present work, the all electron basis set aug-cc-pCV5Z-DK for Be[38] and def2-AQZVPP-JKFI for K[39] are chosen in all calculations. The spectroscopic and molecular parameters are determined with Le Roy’s LEVEL 8.0 program.[40] The PECs calculations are performed with an interval of 0.1 Å over the distance from 2.3 Å to 30 Å for three and two 4Π states and 3.5 Å to 30 Å for three states, and the value decreases to 0.02 Å near Re.

3. Results and discussion
3.1. PECs

The PECs for three , two 4Π, and three electronic states of KBe molecule are performed, which are illustrated in Figs. 13, respectively. We can see that the PECs for all states are smooth. For convenient of discussion, the eight bound electronic states are divided into two categories in the discussion, i.e., the PECs for the doublet () and quartet (4Π, ) electronic states, respectively.

Fig. 1. (color online) PECs for the X2Σ+, 22Σ+, and 32Σ+ states of KBe.
Fig. 2. (color online) PECs for the 14Π and 24Π states of KBe.
Fig. 3. (color online) PECs for the 14Σ+, 24Σ+, and 34Σ+ states of KBe.
3.2. PECs for the states

The dominant electronic configurations at Re for these eight states are tabulated in Table 1. We can see that the dominant compositions for the ground state X2Σ+ at Re is 4σ22228σ2π4; the dominant compositions for the 22Σ+ state can be represented by three main configurations, 4σ22228σ2π4, 4σ22229σ2π4, and 4σ222210σ2π4; the 32Σ+ state can be represented by 4σ22229σ2π4 and 4σ222210σ2π4. The electronic transition and both arise from and electron promotions.

Table 1.

Important electronic configurations of KBe at the equilibrium separations.

.

The PECs for the X2Σ+, 22Σ+, and 32Σ+ states of KBe are plotted in Fig. 1. To our best knowledge, there are not experimental data on spectroscopic parameters of KBe. In previous theoretical studies, Xiao et al.[18] have investigated PECs for the X2Σ+ and 22Σ+ states at the same level, they also given the results of spectroscopic parameters. Compared the results for the X2Σ+ and 22Σ+ states between ours and Xiao’s,[18] the percentage errors are: for the X2Σ+ state, (Ref. [18]), (Ref. [18]), (Ref. [18]), for 22Σ+ state, (Ref. [18]), (Ref. [18]), (Ref. [18]). There are large errors of De and ωe between ours and Xiao’s results, and the reason may be that the 32Σ+ state is considered in the CASSCF and MRCI+Q steps. We also calculate only two states at the same level. Re increased and De decreased for the ground state when the third 3sigma is considered. , , , and The 32Σ+ state may play a very important role, and influence the spectroscopic properties for the X2Σ+ and 22Σ+.

Augustovičvá et al.[9] and Guérout et al.[10] have investigated the ground state for AEM–Alk ( and Sr; , Na, K, Rb, and Cs) at RCCSD(T) and FCI level, respectively, for both MgAlk and SrAlk molecules. Re increases and De decreases with the increase of the atomic weight of Alk, i.e., Re are 3.100, 3.529, 4.087, 4.286, and 4.518 Å; De are 1414.6, 808.4, 637.5, 578.5, and 558.2 cm−1 for LiMg, NaMg, KMg, RbMg, and CsMg,[9] respectively. Re is 3.476, 3.819, and 4.407 Å; De is 2587, 1597, and 1166 cm−1 for LiSr, NaSr, and KSr,[10] respectively. Be, Mg, and Sr in periodic table belongs to group II, the properties of AEM–Alk are similar. Re and De for the ground state of LiBe, NaBe, RbBe, and CsBe molecules are 2.607,[13] 3.052,[13] 3.762,[19] 4.049 Å,[19] and 2340,[13] 1129.2,[13] 761.2,[19] and 418.6 cm−1[19] by Bauschlicher et al. and Yang et al., respectively. The results for the ground state of KBe molecule are as follows: and , respectively. Same to MgAlk and SrAlk molecules, Re(CsBe) > Re(RbBe) > Re(KBe) > Re(NaBe) > Re(LiBe), and De(CsBe) < De(RbBe) < De(KBe) < De(NaBe) < De(LiBe).

The PECs for the 32Σ+ electronic state are studied for the first time. We have given the spectroscopic parameters based on inner electrons correlation. All spectroscopic parameters for the three states are listed in Table 2. Thought there are not experimental data of spectroscopic parameters, we have provided other reference values of spectroscopic parameters for the first three states.

Table 2.

Spectroscopic parameters for eight low–lying states of KBe.

.
3.3. The PECs for quartet states

The dominant electronic configurations at Re for the 14Π and 24Π are 4σ2227σ8σ2π43π and 4σ2227σ9σ2π43π, which is electron promotion for the electronic transition. The 14Π is the ground quartet state, which corresponds to the dissociation channel Be(3P) + K(2P), and has the lowest energy at Re, because the lower dissociation channel Be(1S) + K(2P) can only obtain one doublet Π state.

The dominant electronic configurations at Re for the 14Σ+ are 4σ2227σ8σ9σ2π4 and 4σ2227σ8σ10σ2π4; the dominant electronic configurations at Re for the 24Σ+ are 4σ2227σ8σ9σ2π4, 4σ2227σ8σ10σ2π4, and 4σ2227σ2π43π4π; and the dominant electronic configurations at Re for the 34Σ+ is 4σ2227σ2π43π4π. Therefore, the electronic transition from a state to another at Re can be obtained.

The PECs and spectroscopic parameters for quartet states have been investigated for the first time, the 14Π and 24Π states in Fig. 2 and the 14Σ+, 24Σ+, and 34Σ+ states in Fig. 3. From Fig. 3, it can be seen that the avoided crossing of 24Σ+ and 34Σ+ states occurred at about 6.0 Å. The energy difference between these two states is only 1.909 cm−1 at . All spectroscopic parameters for two 4Π and three states are listed in Table 2. De for the 14Σ+, 24Σ+, and 34Σ+ states are 24.927, 98.330, and 22.178 cm−1, respectively. They are shallow potential wells.

3.4. Molecular parameters

Based on the reliable PECs for the bound states (X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+), we have determined their vibrational levels Gυand inertial rotation constants Bυ when . Gυ and Bυ for the ground state together with Xiao’s results[18] are collected in Table 3, the other data are collected in Supplementary Material. Compared with the results of Xiao et al.,[18] we overestimate the Gυ for the ground state. From Supplementary Material, the first 21 levels of Gυ for the 22Σ+ state are close to the results of Xiao et al.[18] The largest error is only about 2.6%, which is in agreement with the results of De.

Table 3.

The vibrational levels and inertial rotation constants Bv for the ground state X2Σ+ of KBe.

.

The highest vibrational level for the ground state X2Σ+ is in our work. This level lies 0.7035 cm−1 below the dissociation limit, or the energy of 99.91% for De. The level is for the 22Σ+ state, which lies 0.013 cm−1 below the dissociation limit, the energy of almost 100% for De. The highest vibrational level for the ground state X2Σ+ obtained by Xiao et al.[18] is , which lies 1.02 cm−1 below the dissociation limit, the energy of 99.85% for De.

Inertial rotation constants Bυ for states (X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+) are also listed in Supplementary Material. Based on the accuracy of calculations, we believe that the values of Gυ and Bυ for these eight states are reliable.

3.5. Dipole moments, transition dipole moments, and radiative lifetimes

The dipole moment as a function of R is plotted in Fig. 4, for the X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+ states. It is clearly seen that the avoided crossing position between the 24Σ+ and 34Σ+ states occurs at about 6.0 Å. The dipole moments for eight states approach zero at large interatomic distance. We can say that the dissociation channels for these eight states are the sum of neutral K atom and Be atom. Dipole moments at Re for eight states are listed in Table 4.

Fig. 4. (color online) Dipole moments for the X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+ states of KBe.
Table 4.

Dipole moments at the equilibrium distance for the eight states of KBe.

.

Notable, for one-electron and two-electron systems, the dipole moments for the ground state at Re for KBe and LiYb[22] are 1.748 and 0.03 Debye (D), respectively; and the electronegativities (EN) of Li, K, Be, and Yb atoms are 0.98, 0.82, 1.57, and 1.06 eV, respectively. It is clearly seen that the difference of EN between K and Be is 0.75 eV, but the value is only 0.08 eV between Li and Yb. Therefore, KBe has a sizable dipole moment, but LiYb does not. For Alk-Be series, the dipole moments are 2.415,[14] 1.978,[15] 1.748, 1.302,[19] and 1.046 D[19] for the ground state of LiBe, NaBe, KBe, RbBe, and CsBe molecules, respectively. It can be seen that the dipole moments for the ground state decrease with the increase of the atomic weight of Alk.

Figure 5 shows the transition dipole moment functions (TDMs) for the 22Σ+–X2Σ+, 32Σ+–X2Σ+, and 24Π−14Π systems of KBe molecule. Notable the transitions for the 22Σ+–X2Σ+ and 24Π−14Π+ have the same asymptotic limit, and come from an allowed transition of K atom at large interatomic distance. Our calculation results at are 7.8185 and 7.8124 D, respectively. Xiao et al.[18] obtains the value 8.17 D for the 22Σ+–X2Σ+ bond system, which is about 0.25 D bigger than our results. The 32Σ+––X2Σ+ system is transition of Be atom at large interatomic distance, which is a forbidden transition. The TDMs for the 32Σ+–X2Σ+ system decreases and tends to be zero at large interatomic distance.

Fig. 5. (color online) Transition dipole moments for the 22Σ+–X2Σ+, 32Σ+–X2Σ+, and 24Π−14Π transitions of KBe.

Transition probabilities for emission including Einstein coefficients Aυ′υ″ and Franck–Condon factors qυ′υ″ for the 22Σ+–X2Σ+, 32Σ+–X2Σ+, 32Σ+−22Σ+, and 24Π−14Π band systems have been evaluated, which are tabulated in Table 5. The values for transitions from the vibrational levels of a higher state to the vibrational levels of a lower state are listed in Table 5. For the first five vibrational levels () for the 22Σ+–X2Σ+ system, the stronger transitions occur from to , 1, 2, 3; to , 2, 3, 4, 5, 6; to , 1, 2, 6, 7, 8; to , 1, 4; and to , 3. We can clearly see that both Aυ′υ″ and qυ′υ″ are concordant on this prediction for the calculated transition systems.

Table 5.

Einstein spontaneous emission coefficients Aυ′υ″ (in s−1) and Franck–Condon factors qυ′υ″ (in italics) for the , , , , and band systems of KBe.

.

The radiative lifetimes can be calculated as the reciprocal of total Einstein spontaneous emission coefficient,

The lifetime of the first 11 vibrational levels for the 22Σ+–X2Σ+, 32Σ+–X2Σ+, and 24Π−14Π transitions of KBe are listed in Table 6. We can see that three transitions are strong. When , the lifetimes are 61.28, 87.01, and 71.03 ns for the 22Σ+–X2Σ+, 32Σ+–X2Σ+, and 24Π−14Π transitions, respectively.

Table 6.

Radiative lifetime (in ns) of the transitions 22Σ+–X2Σ+, 32Σ+–X2Σ+, and 24Π−14Π band systems of KBe molecule.

.
4. Conclusion and perspectives

In this paper, the potential energy curves (PECs) for electronic states (X2Σ+, 22Σ+, 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+) of KBe molecule are performed by using highly accurate multireference configuration interaction (MRCI) with aug-cc-pCV5Z-DK for Be atom and def2-AQZVPP-JKFI for K atom.

The spectroscopic parameters for these states are also predicted. The 14Π is the ground quartet state, three quartet Σ+ electronic states are shallow potential wells. Based on the accurate results of PECs, molecular parameters for these eight states are also obtained by solving the radial Schrödinger equation. The dipole moments, transition dipole moments, Franck–Condon factors qυ′υ″, Einstein coefficients Aυ′υ″ and the radiative lifetimes τυ′ for the 22Σ+–X2Σ+, 32Σ+–X2Σ+, and 24Π−14Π transitions are also predicted at the same time. The 32Σ+, 14Π, 24Π, 14Σ+, 24Σ+, and 34Σ+ states are investigated for the first time.

Reference
[1] Jones K M Tiesinga E Lett P D Julienne P S 2006 Rev. Mod. Phys. 78 483
[2] Wang D Qi J Stone M F Nikolayeva O Wang H Hattaway B Gensemer S D Gould P L Eyler E E Stwalley W C 2004 Phys. Rev. Lett. 93 243005
[3] Takasu Y Komori K Honda K Kumakura M Yabuzaki T Takahashi Y 2004 Phys. Rev. Lett. 93 123202
[4] Kitagawa M Enomoto K Kasa K Takahashi Y Ciurylo R Naidon P Julienne P S 2008 Phys. Rev. A 77 012719
[5] Nemitz N Baumer F Münchow F Tassy S Görlitz A 2009 Phys. Rev. A 79 061403
[6] Kerman A J Sage J M Sanis S Bergeman T Demille D 2004 Phys. Rev. Lett. 92 033004
[7] Kerman A J Sage J M Sanis S Bergeman T Demille D 2004 Phys. Rev. Lett. 92 153001
[8] Lane I C 2012 Phys. Chem. Chem. Phys. 14 15078
[9] Augustovičová L Soldán P 2012 J. Chem. Phys. 136 084311
[10] Guérout R Aymar M Dulieu O 2010 Phys. Rev. A 82 042508
[11] Gopakumar G Abe M Kajita M Hada M 2011 Phys. Rev. A 84 062514
[12] Chen K Yang C L Wang M S Ma X G Liu W W 2012 Spectrochim. Acta A 99 57
[13] Bauschlicher C.W. Jr. Langhoff S R Partridge H 1992 J. Chem. Phys. 96 1240
[14] Svetlana K Alexander P Maria L Jacek K Paul S J 2011 J. Chem. Phys. 135 164108
[15] Pablo J B Friedrich G 2002 Mol. Phys. 100 1681
[16] Xiao K L Yang C L Wang M S Ma X G Liu W W 2013 J. Quant. Spectrosc. Radiat. Transf. 129 8
[17] Bruna P J Grein F 2002 Mol. Phys. 100 1681
[18] Xiao K L Yang C L Wang M S Ma X G Liu W W 2013 J. Chem. Phys. 139 074305
[19] You Y Yang C L Wang M S Ma X G Liu W W Wang L Z 2015 J. Quant. Spectrosc. Radiat. Transf. 165 56
[20] Tohme S N Korek M 2013 Chem. Phys. 410 37
[21] Gopakumar G Abe M Das B P Hada M Hirao K 2010 J. Chem. Phys. 133 124317
[22] Zhang P Sadeghpour H R Dalgarno A 2010 J. Chem. Phys. 133 044306
[23] Wan M J Shao J X Gao Y F Huang D H Yang J S Cao Q L Jin C G Wang F H 2015 J. Chem. Phys. 143 024302
[24] Wan M J Shao J X Huang D H Jin C G Yu Y Wang F H 2015 Phys. Chem. Chem. Phys. 17 26731
[25] Wan M J Huang D H Shao J X Yu Y Li S Li Y Y 2015 J. Chem. Phys. 143 164312
[26] Werner H J Knowles P J 1988 J. Chem. Phys. 89 5803
[27] Knowles P J Werner H J 1988 Chem. Phys. Lett. 145 514
[28] Li R Wei C L Sun Q X Sun E P Jin M X Xu H F Yan B 2013 Chin. Phys. B 22 123103
[29] Zhu Z L Qiao H Lang J H Sun J F 2013 Chin. Phys. B 22 103102
[30] Wang J M Feng H Q Sun J F Shi D H 2012 Chin. Phys. B 21 023102
[31] Liu X J Miao F J Li R Zha C H Li Q N Yan B 2015 Acta Phys. Sin. 64 123101 in Chinese
[32] Xing W Liu H Shi D H Sun J F Zhu Z L 2016 Acta Phys. Sin. 65 033102 in Chinese
[33] Werner H J, Knowles P J, Lindh R, et al. MOLPRO, version 2010.1 http://www.molpro.net
[34] Knowles P J Werner H J 1985 J. Chem. Phys. 82 5053
[35] Knowles P J Werner H J 1985 Chem. Phys. Lett. 115 259
[36] Werner H J Knowles P J 1988 J. Chem. Phys. 89 5803
[37] Knowles P J Werner H J 1988 Chem. Phys. Lett. 145 514
[38] Prascher B Woon D E Peterson K A Dunning T H Jr. Wilson A K 2011 Theor. Chem. Acc. 128 69
[39] Weigend F 2008 J. Comput. Chem. 29 167
[40] Roy R J Le. LEVEL 8.0: A Computer Program for Solving the Radial Schrödinger Equation for Bound and Quasibound Levels Chemical Physics Research Report CP-663 University of Waterloo 2007